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ABSTRACT 

We propose a plasma experiment to be used to investigate fundamental properties of astro- 
physical dynamos. The highly conducting, fast-flowing plasma will allow experimenters to explore 
systems with magnetic Reynolds numbers an order of magnitude larger than those accessible with 
liquid-metal experiments. The plasma is confined using a ring-cusp strategy and subject to a 
toroidal differentially rotating outer boundary condition. As proof of principle, we present mag- 
netohydrodynamic simulations of the proposed experiment. When a von Karman-type boundary 
condition is specified, and the magnetic Reynolds number is large enough, dynamo action is 
observed. At different values of the magnetic Prandtl and Reynolds numbers the simulations 
demonstrate either laminar or turbulent dynamo action. 

Subject headings: MHD, plasmas, stars: magnetic fields, ISM: magnetic fields 



1. Introduction 

One of the central problems of astrophysi- 
cal fluid dynamics is the magnetohydrodynamic 
(MHD) dynamo, the process by which the mo- 
tion of an electrically conducting fluid amplifies a 
small seed magnetic field until the field becomes 
dynamically important. Astrophysical dynamos 
can be categorized into two types: the large-scale 
dynamo, where the scale of the generated mag- 
netic field is the same as or larger than the system 
generating it, such as in planets or stars, or the 
small-scale dynamo, where the scale of the mag- 
netic field is much smaller than the sys tem size, 
such as in warm interstellar medium (jKulsrud 
1999). Of particular physical importance in this 
context is the system's magnetic Prandtl number, 



Pm = v/rj, where v is the kinematic viscosity 
of the fluid and r\ its magnetic diffusivity, whose 
value indicates the scale of the velocity field at 
the resistive-diffusive cutoff. When Pm -C 1 the 
magnetic-diffusion scale is in the hydrodynamic 
inertial range and the primary interactions are at 
that scale. If a dynamo occurs in this situation 
it is often large scale, though recent evidence in- 
dicates t hat a small-scale dynamo may also be 



possi ble (jlskakov et al.l 120071 : iSchekochihin et al 



2007). Conversely, Pm 3> 1 gives a magnetic- 



diffusion scale below the viscous-diffusion cutoff, 
permitting the stretching and folding of the mag- 
netic field at small scales and allowing the oppor- 
tunity for a small-scale dynamo. 

The small-scale dynamo has only been studied 
theoretically and through numerical simulations 
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see lHaugen et al.l ([20041 ): ISchekochihin et al 
(J2007I) and references therein). The simulations 
tend to be done in periodic boxes with random 
non-helical forcing, in the absence of a mean 
flow, allowing the study of the dynamo under ho- 
mogeneous and isotropic fluctuating conditions. 
Most simulations have been done in the range of 
Pm ~ 1, meaning with identical viscous and mag- 
netic dissipation scales, though there have been 
not a ble exceptions where much higher values of 



mghc r values o: 
Y3, iKinnev et al 



Pm (ISchekochihin et al.ll2004l200 

2000) and values with Pm < 1 (Isk akov et al 



20071 : ISchekochihin et all l2007h have been used. 
Despite the large body of numerical and theo- 
retical work on such dynamos, small-scale dy- 
namos have never been experimentally studied 
due to the dearth of fluids for which Pm ^> 1 or 
Pm ~ 1 . The only fluids which satisfy these crite- 
ria are plasmas, and plasma dynamo experiments, 
to study large or small-scale dynamos, have not 
thus far been constructed. 

In contrast, several groups have recently been 
successful in achieving large-scale dynamo action 
in experiments using liquid sodium, a fluid for 
whic h Pm«l dGailitis et al.ll2000l :l stiegritz fc Miiller 



20011 iMonchaux et al.ll2007l ). In such experiments 
the magnetic Reynolds number, Rm — v^a/m 
where vq is a characteristic speed, and a a length 
scale, characterizes the ratio of magnetic field ad- 
vection to diffusion. For idealized laminar veloc- 
ity fields the critical value of Rm for magnetic 
self-e xcitation, Rm rr \t, i s predicted to be a round 
100 (|Forest et al.l l2002t iRavelet et al.l l2005h . To 
achieve this using liquid sodium in an experiment 
of radius a = 0.5 m requires a mechanical input 
power of Pmech ~ 100 kW. These flows are very 
turbulent, however, and turbulent flows have the 
unfavorable scaling Rm cx (Pmechfl) 1 ^ 3 - It would 
be interesting to study dynamo physics at a higher 
order of magnitude, Rm — 1000 for example, but 
to achieve such a value of Rm in a sodium exper- 
iment of similar size would require a mechanical 
input power of 100 MW! This is a serious limi- 
tation to addressing higher- Rm dynamo regimes 
using liquid sodium. 

A different class of fluids which could be used 
to study dynamo physics is plasma. To match 
the conductivity of sodium a singly ionized plasma 
requires an electron temperature of 630 eV, a 
plasma temperature only found in fusion exper- 



iments. However, plasma flows can be efficiently 
driven to much higher speeds than liquid metals. 
Thus for a plasma experiment of similar size (a = 1 
m) to achieve a value of Rm = 1000 would simul- 
taneously require an electron temperature of only 
T c = 10 eV and a velocity vq = 20 km s -1 . These 
are modest values which can easily be achieved in 
many plasma confinement configurations. Plasma 
experimenters also have the novel ability to change 
the viscosity of their plasmas by several orders of 
magnitude. Such flexibility could allow the ex- 
ploration of laminar and turbulent dynamos in a 
wide range of Pm, from Pm < 1 to Pm 1. 
Such variability would be a striking improvement 
over liquid- metal experiments. 

However, plasmas suffer from an important dis- 
advantage when compared to liquid metals: they 
can be difficult to control and confine without the 
presence of a magnetic field. This is a problem for 
the study of dynamos, of course, because the pres- 
ence of a dynamically significant magnetic field 
changes the fundamental behavior of the problem 
under study. Those challenges common to most 
plasma experiments-the need for thermal confine- 
ment to keep the plasma away from container walls 
and hot enough to be a good conductor, and some 
scheme to drive the flow-must be overcome, in 
a dynamo experiment, without magnetizing the 
plasma. A carefully designed confinement scheme 
is needed to meet such criteria. 

In this paper we introduce a multipurpose 
plasma experiment propo sed by one of the au- 



thors (jForest et al.l l2008h . The experiment is 



spherical and based on an axisymmetric ring- 
cusp confinement strategy. The confining mag- 
netic field is localized to the periphery of the ex- 
periment and a large, unmagnetized plasma vol- 
ume is created in the experiment's core. The 
only considered means of injecting energy into 
the plasma, once formed, is by forcing the velocity 
field in the toroidal direction at the plasma's outer 
edge. Unlike spherical Couette flow experiments, 
whose outer boundary rotates as a rigid rotor, the 
proposed experiment's outer boundary is capa- 
ble of differential forcing as a function of angle, 
meaning, approximately, a velocity field subject to 
the boundary condition v(r = a,9,<p) — v<f,(6)<p, 
where a is the radius of the sphere. The only 
physical restriction on the boundary condition 
is that it remain zero at the poles, meaning 
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u^(0) = v^(%) — 0. By adjusting the plasma 
such that Pm = 1.0, and choosing Rm appropri- 
ately, the experiment should be able to generate 
a laminar dynamo, conditions that have never be- 
fore been created in the laboratory. It should also 
be able to produce turbulent dynamos by oper- 
ating high-fluid- Reynolds-number flows with high 
values of Rm. It would also not be restricted to 
dynamo physics, as a vast set of velocity fields 
could be generated in a simply-connected volume, 
allowing the experimental study of numerous as- 
trophysical and geophysical phenomena. Further 
details of the proposed experiment are presented 
in Section [2] 

We also present numerical simulations as proof- 
of-concept of the proposed experiment, simulating 
a few of the many possible configurations us- 
ing a single-fluid approximation of the plasma. 
The simulations are performed using a paral- 
lelized version of a three-dimensional incompress- 
ible non-linear MHD code develope d to simulate 



the Madison Dynamo E xperiment (jBavliss et al 



20071 : iReuter et ai]l2008l ). By varying the outer 



toroidal velocity field boundary condition differ- 
ent flow regimes have been studied. In Section [3] 
we present a boundary condition which results 
in flows which display laminar or turbulent dy- 
namo action depending on the simulation's value 
of the magnetic Reynolds and Prandtl numbers. 
We conclude with a discussion of possible avenues 
of research for the experiment. 

2. Experimental description 

The plasma is confined in a spherical ge- 
ometry using an axisymmetric ring-cusp strat- 
egy, consisting of rings of permanent magnets 
mounted to the inner wall of the sphere. The 
poles of the magnets are oriented radially, and 
the rings have alternating polarity. A partial 
cross section of this configuration is presented 
in Figure [1] Such a confinement strategy has 
bee n used previously in a cylindrical geome 
try ( Limpacc her fc MacKenzie _ 1973 1 Leung et al 



197alLang fc Hershkowit Jl978tlCho et alll98gj k 
and is currently used in ion sources f or neu 
tral beam heating ( Ehlers fc Leune 1979h and in 



the c ontext of plasma processing (jPelletier et al 



within a radial distance on the order of the space 
between the magnets. The experimental volume 
is essentially magnetic field free, resulting in a 
very high-/3 plasma, i.e. pv 2 3> B 2 /po, where p 
is the plasma mass density and po the magnetic 
permeability of a vacuum. 

Confinement in such a magnetic geometry has 
been well studied. Particle confinement is deter- 
mined by ion-acoustic flows into the cusps, which 
can be modeled as a loss area equal to the lin- 
ear d imension of the cusp tim es the ion gyrora- 
dius (jHershkowitz et al. 1 119751 ). Electrons are well 



confined by the cusp fields ( Leung et al. 1976h : 



their primary energy losses are through the cusp 
and through ion collisions. Ion energy confinement 
is mostly determined by charge exchange losses re- 
sulting from the relatively high fraction of neutrals 
in such devices. The ionization fraction of previ- 
ous and current multidipole confinement systems 
varies considerably, but is typically < 0.2. 

Improved technology could allow the proposed 
experiment to exceed the confinement limits of 
earlier apparatuses. Modern NeFeB permanent 
magnets can generate cusp fields twice as strong 
as those in previous cylindrical experiments. This 
would result in lower loss rates, due to the narrow- 
ing of the cusps, and increased electron temper- 
ature, which has been sh own to scale with cus p 
magnetic field strength (jEhlers fc Leune 1982). 
Heating of the plasma could be done using large 
area Lanthanum Hexaboride (LaBg) cathodes. 
Such ca thodes have previously been used in ion 
sources dEhlers fc Leundll979i iLeung et alJll984t 
Fmcosv fc Leundll98"flK and in a variety of plasma 



19841 ). The magnetic field is localized to the outer 
edge of the sphere and drops to a negligible value 



experiments ( Ono et aDll987 : Darrow et alJ 1990). 
These sources produce significantly more power 
than tungsten filaments, which would lead to 
higher temperatures and ionization fractions than 
previous devices. 

To inject momentum into the experiment, ring 
anodes and cathodes are alternately placed be- 
tween the ring magnets (Figure [1]) . These gen- 
erate an electric field that causes an E X B drift 
in the toroidal direction. The magnitude and di- 
rection of this drift can be varied as a function of 
angle by changing the electrical potential between 
the anodes and cathodes, allowing control over the 
fluid's outer toroidal velocity field boundary con- 
dition. Viscous coupling between the magnetized 
and unmagnetized regions of the experiment is as- 
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Fig. 1. — Partial vertical cross section of the pro- 
posed experiment. Rings of permanent magnets, 
of alternating polarity, line the inside of the sphere 
with their poles oriented radially. Ring anodes and 
cathodes lie between the magnets. The resulting 
E X B drift is in the toroidal direction. By varying 
the potential between the anodes and cathodes the 
forcing at the outer boundary can be customized. 
The apparatus is axisymmetric. 



sumed. 

A broad set of experimental regimes can be gen- 
erated using such a device. Some of the specifi- 
cations of the proposed experiment are listed in 
Table Q] In the context of this study, we are par- 
ticularly interested in the magnetic Reynolds and 
Prandtl numbers. The magnetic Reynolds number 
for a plasma is given by 

J2m= 1.5^5^, (1) 
Z 

where the electron temperature, T , is measured in 
electron volts, U is the peak speed of the plasma 
in km s , the length scale a is measured in me- 
ters, and Z is the charge of the ions. High electron 
temperatures, and the high speeds of plasmas, can 
lead to very large values of magnetic Reynolds 
number. In an unmagnetized plasma the viscosity 
is given by v ~ VyjT e j, where Vt\ — 2ylimi is the 
thermal velocity of the ions, t g ; is the electron- 
ion collision time, T; is the ion temperature and 
mi is the ion mass. The unmagnetized magnetic 
Prandtl number is then 

j,3/2rp5/2 

Pm = 0.18 c „ 1 — , (2) 

where the ion temperature is measured in electron 
volts, \x is the atomic mass number of the ions, and 
n is the number density in units of 10 18 m~ 3 . The 
strong dependence on density and ion mass allows 
a very large range of magnetic Prandtl number, 
and by extension fluid Reynolds number, Re = 
Rm/Pm, allowing the experimenter the ability to 
specify whether a given plasma will be laminar or 
turbulent. 

Varying density and ion mass to control the vis- 
cosity of the plasma is not without its tradeoffs, 
the primary one being the saturation level of the 
magnetic field generated by the dynamo, since it 
is dependant upon the kinetic energy of the fluid. 
Assuming that saturation occurs when the mag- 
netic energy is in equipartition with the kinetic 
energy, B 2 /Mo = pU 2 , implies 

Bequipartion ~ QA^fjmU GauSS. (3) 

At T e = 10 eV, U = 20 km s" 1 , ^cquipartion could 
be as high as 240 Gauss with Ar at n = 10 19 
m . However, as we will show in Section [3] 
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Table 1 

Parameters of the Experiment 



Quantity 


Symbol 


Value 


Unit 


plasma radius 


a 


1.5 


m 


number density 


n 


10 17 -10 19 


m -3 


ion temperature 


Ti 


0.5-4 


eV 


electron temperature 


T c 


2-10 


eV 


peak speed 




0-20 


km s _1 


ion species 


H, He, Ar 


1, 4, 40 


amu 


pulse length 


Tpulsc 


5 


s 


plasma beta 


P 


10 4 




resistive time 

magnetic Reynolds number 
Reynolds number 
magnetic Prandtl number 


JL?Ti max 

Re 
Pm 


50 

- 1000-2000 
2.4xl0 1 -3.8xl0 6 
S.OxlO^^.GxlO 1 


ms 



simulations indicate that in saturation the mag- 
netic field magnitude is not this high, but at best 
B 2 /fi - O.lpU 2 . 

The lack of magnetic field in the volume of the 
experiment and a method to control the velocity 
field means that this device should satisfy the cri- 
teria needed to create magnetically self-exciting 
plasmas. Its high value of Rm and variable mag- 
netic Prandtl number would allow it to be used to 
study physical regimes inaccessible to liquid-metal 
experiments. It would be the first plasma experi- 
ment used to stu dy dynamos thou gh not the first 
to be proposed ( Wang et al.l 12002^) . and its wide 
range of parameters (Table [T]) should allow it to 
be useful for the study of a number of other as- 
trophysical phenomena: magnetorotational insta- 
bility, high-/3 instabilities and rotating convection, 
to name just a few. In the next section we present 
numerical simulations which demonstrate that this 
concept should succeed as a dynamo experiment. 
The other topics above will be examined in future 
work. 

3. Numerical simulations 

A comprehensive simulation of such a plasma 
experiment would consider both the electron and 
ion species, as well as include a careful treat- 
ment of the E X B forcing at the outer bound- 
ary. As a first approximation, however, we treat 
the plasma as an incompressible single fluid. This 



is a justified approximation since the mean free 
path of the particles is much shorter than the sys- 
tem size, and sound waves are not relevant. The 
simulation solves the vorticity evolution and mag- 
netic induction equations using a standard pseu- 
dospectral method based on spherical harmon- 
ics, truncated at maximum spherical harmonic de- 
grees t 



m ax and m ma ,x; details have been given pre - 



viously (|Bavliss et all 120071 : iReuter et alj 12008) . 



The E X B forcing is approximated as a toroidal 
velocity field non-zero no-slip boundary condition. 
An electrically-insulating outer boundary is as- 
sumed. 

For the sake of brevity, here we restrict our 
attention to two values of magnetic Reynolds 
number, Rm — 300, 2000, where the magnetic 
Reynolds number is based on the peak speed of the 
toroidal outer boundary condition, which is unity 
in all simulations. We examine several values of 
Pm to demonstrate that this experiment should 
be able to generate both laminar and turbulent dy- 
namos. Future studies will examine other regimes 
of interest. 

3.1. Toroidal boundary conditions 

To demonstrate the flexibility of such an ap- 
paratus, four different axisymmetric steady-state 
velocity fields, generated using different outer 
boundary conditions with Rm = 300 and Pm = 
1.0, are presented in Figure [2] The first, perhaps 
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the simplest that can be generated in such a de- 
vice, is a rigid-rotor-type azimuthal field. Such a 
velocity field can be used, in combination with the 
injection of light ions into a heavy plasma, to study 
rotating convection in plasmas. The second veloc- 
ity field, Figure [2jb), is generated by a boundary 
condition that follows a ~ s _1 dependence for 
much of its angular extent, where s is the cylindri- 
cal radial coordinate. The toroidal flow generated 
by this boundary does not have a s _1 dependence, 
but rather is Keplerian for much of the volume of 
the sphere, with ~ s -1 / 2 at the equator, and 
as such should be unsta ble to the magnetorota - 
tional instability (MRI) (jBalbus fc Hawlevlll99ih . 
Not surprisingly, when this flow is exposed to an 
externally-applied magnetic field an MRI-like in- 
stability appears to develop. 

The velocity fields presented in Figures [2jc) 
and [2jd) are both magnetically unstable, dis- 
playing dynamo action under the conditions pre- 
sented. The first boundary condition is equato- 
rially symmetric, resulting in a flow with net ki- 
netic helicity, while the second boundary condi- 
tion is equatorially anti-symmetric. For the re- 
mainder of this paper we will focus on the sec- 
ond of these t wo cases, which is b ased on the von 
Karman flow (jvon Karmanl Il92lh , a flow which 
has been extensively studied in the cylindrical ge- 
omet ry (|Odier et alJ Il998t iBourgoin et"aD I2002L 
20041 ). In this case the outer boundary rotates 
in opposite directions near the poles of the sphere 
and relatively little near the equator, as seen in 
Figure [3l The boundary condition is constructed 
from non-zero outer boundary values for the even- 
numbered axisymmetric spherical harmonic com- 
ponents I = 2,4,6,8; these values, and all the 
boundary values used in this study, are given in 
Appendix [X] This boundary condition results in 
laminar and turbulent dynamo action depending 
on the value of Rm and Re (for this boundary 
condition the transition to turbulence begins at 
Re ~ 500, depending on the level of forcing). The 
parameter values used in the simulations presented 
in this study, and exemplary physical values of cor- 
responding plasmas, can be found in Table [2] 




Fig. 2. — Steady-state axisymmetric velocity fields 
generated by different toroidal boundary condi- 
tions, with Rm = 300 and Pm = 1.0. The axis of 
symmetry runs vertically. In the left hemisphere 
are the contours of toroidal speed, and in the right 
hemisphere are the contours of the poloidal stream 
function. The boundary conditions are (a) solid- 
body rotation, which generates no poloidal flow 



(b) 



which results in a Keplerian profile 



(c) equatorially symmetric, which results in dy- 
namo action and (d) von Karman (given in Fig- 
ure^, which also generates a dynamo. Note that, 
as indicated in Figure El the peak speed is set to 
1.0, but the scale range in figures (b)-(d) has been 
reduced for clarity. 



3.2. Laminar dynamo simulations 

Among the many appealing features of the pro- 
posed experiment is its ability to create laminar 
MHD flows at very low values of fluid Reynolds 
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Table 2 

Simulation Parameters & Plasma Values 









T c 


Ti 


U 


n 












Rm 


Pm 


Re 


(eV) 


(eV) 


(km s- 1 ) 


(10 18 m- 3 ) 


Ion 


(G) 


N r 


'-max 


^ m ax 


300 


2.0 


150 


4.3 


3.1 


20.0 


13.8 


H 


3.2 


200 


14 


14 


300 


1.0 


300 


4.3 


3.1 


20.0 


27.6 


H 


4.2 


200 


14 


14 


300 


0.5 


600 


10.0 


4.0 


5.7 


22.7 


He 




300 


18 


14 


2000 


1.0 


2000 


15.9 


4.0 


19.1 


23.0 


He 


12.0 


600 


50 


30 



Note. — The radius of the unmagnetized plasma volume, a, is taken to be 1.1 m. N r is the 
number of radial grid points, equally spaced in the range < r < 1. B max refers to the maximum 
saturated magnetic field strength. 



to 



> 




1.0 1.5 2.0 

Angle [Radians] 

Fig. 3. — Toroidal boundary condition used in the 
simulations to generate a von Karman-type flow. 
This boundary condition is constructed from the 
boundary values specified in Appendix [A"l 



number, something that is not possible with 
liquid- metal experiments. Though such flows are 
an imperfect analogy of naturally-occurring dy- 
namos, since all natural flows are turbulent, their 
study can be used to test basic MHD under ide- 
alized conditions as well as the numerical codes 
which simulate them. Low- Reynolds-number flows 
can also be used to study the laminar-to-turbulent 
transition, and should allow the observation of 
the spontaneous relaminarization o f a turbulent 



flow by a self-excited magnetic field (jBavliss et al 



2007). Here we restrict our examination to a set 



of laminar dynamos which could be generated in 
the proposed experiment. 

It might come as a bit of a surprise that a 
laminar velocity field generated by a toroidal 
differentially-rotating outer boundary can gen- 
erate a dynamo, since the poloidal component of 
the velocity field only results from Ekman circu- 
lation, making it difficult to achieve the balance 
of toroidal-to-poloidal flow needed for laminar dy- 
namo action. Also, since the toroidal flow peaks at 
the edge of the sphere, the poloidal field is unable 
to easily stretch the magnetic field around the 
toroidal flow's p eak, and thus the stretch- twist- 
fold mechanism (jChildress fc Gilbertl I1995I ) that 
sustains laminar dynam os is much less efficien t 
than that of other flows ( Dudley k. Jamesl 1989h . 
This manifests itself in the comparably large crit- 
ical magnetic Reynolds number required for the 
toroidally-driven flows to self-excite. 

The velocity field which results from the von 
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Karman boundary condition, with Rm — 300 
and Pm = 1.0 (Figure H[d)), is axisymmetric, 
counter-rotating in the toroidal direction, and has 
a poloidal flow which rolls inward at the equator 
and outward at the poles. As such the velocity 
field is qualitatively si milar to the s2t2 flow of 
Dudley k. Jamed ( 19891 ). with the important dis- 
tinction that the toroidal flow peaks at the sphere 
boundary. The magnetic and kinetic energies of 
this simulation, as a function of time in resis- 
tive units (t v — a 2 ft]), are given in Figure [4] 
The growth rate of the magnetic energy is very 
small, symptomatic of the flow's inefficiency as 
a dynamo. The critical magnetic Reynolds num- 
ber for this flow, based on a linear stability anal- 
ysis, is i?m cr it ~ 237. As is required for ax- 
i symmetric vel ocity fields by Cowling's theorem 
(ICowlindll93"3h . 1-he excited magnetic field is non- 
axisymmetric, dominated by m = 1 modes; this is 
a large-scale dynamo. Overall the final magnetic 
field is weak, with the magnetic energy peaking 
an order of magnitude below the kinetic energy, 
resulting in minimal modification of the velocity 
field during saturation. The large-scale structure 
of the saturated magnetic field is given in Fig- 
ure O demonstrating a geometry si milar to other 
s2£2-type dynamo magneti c fields (iBavliss et aT 



20071 : iGissinger et~al]l2008al |bl: iReuter et ail 2009) . 



Using the parameters in Table El the value of 
the magnetic field contour in Figure [5] is about 
2 Gauss, which should be measurable in the pro- 
posed experiment using Hall-effect sensors. 

The experimenter's ability to vary the fluid's 
magnetic Prandtl number can be used to explore 
the onset of turbulence and how that turbulence 
affects laminar dynamo action. The energies of 
two other simulations which use the boundary 
condition presented in Figure [3l with Rm — 
300 but with Pm = 2.0,0.5 (Re = 150,600), 
are plotted in Figure [U For Pm = 2.0, the 
boundary condition results in a velocity field 
that is magnetically unstable. This flow has a 
higher critical magnetic Reynolds number than 
the Pm = 1.0 case, Rm cri t — 274, and conse- 
quently has a lower magnetic field growth rate, 
since as expected the growth rates are propor- 
tional to Rm — Rm cr i t . Interestingly, the mag- 
netic field energy saturates at essentially the same 
value as the Pm = 1.0 case, indicating that 
the saturated magnetic field does not follow the 



CO 



10° 
10" 5 
10" 10 
10" 15 
10' 20 
10" 25 

1Q -30 



V 

B, P m =0.5 
B, P m =1.0 
B, P m =2.0 







10 



15 20 
Time [<cj 



25 



30 35 



Fig. 4. — Kinetic and magnetic energies versus 
time for the system forced by the boundary con- 
dition given in Figure El with Rm = 300 and 
Pm = 0.5,1.0,2.0. The kinetic energy for the 
three simulations is essentially the same. The 
Pm = 0.5 simulation does not display dynamo 
action. 
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Fig. 5. — Contour of constant magnetic field mag- 
nitude during saturation, for the Rm — 300, 
Pm =1.0 laminar dynamo. The axis of symmetry 
runs vertically. 
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B 2 ~ [pv I ! aa 2 )(Rm — Rrn rr \ f ) scaling, as migh t 
otherwise be expected (jPetreris &: Fauvd l200lh . 
The Pm = 0.5 case does not magnetically self- 
excite, with the energy used to initialize the mag- 
netic field modes quickly dissipating away. In this 
case the flow contains several non-axisymmetric, 
though stationary, components which interfere 
with the growth of the magnetic field. This 
underlying symmetry-breaking hydrodynamic in- 
stability is a first step towards developed tur- 
bulence, and as such permits the study of how 
non-axisymmetric contributions change the onset 
conditions for dynamo action. 

3.3. Turbulent dynamo simulations 

The question of the role of turbulence in the 
development and maintenance, or hindrance and 
destruction, of astrophysical magnetic fields is of 
ongoing importance. This is especially true in 
light of the recent tur bulent-dynamo results from 
the VKS2 experiment ( Monchaux et al.ll2007 ) and 
the suggestion that coherent turbulence may be re - 
sponsible for that dynamo ( Laguerre et al. 2008h . 
Clearly the proposed experiment will be most rel- 
evant to the study of astrophysical dynamos if it is 
able to generate self-excited magnetic fields under 
turbulent conditions. 

To examine this question, simulations were per- 
formed at the higher end of the experiment's ex- 
pected range of magnetic Reynolds number, Rm = 
2000, with Pm = 1.0 (Re = 2000), using the 
toroidal outer boundary condition presented in 
Figure [3l The kinetic and magnetic energies of 
this simulation versus time are presented in Fig- 
ureJHl The magnetic energy grows exponentially in 
time, saturating an order of magnitude below the 
kinetic energy. Using the parameters in Table [2] 
the magnetic field in saturation peaks as high as 
12 Gauss. Since the magnetic field varies rapidly 
in time it should be easily measurable using B-dot 
coils. 

Both the velocity and magnetic fields fluctuate 
wildly throughout the simulation. Snapshots of 
contours of |B| during saturation are presented in 
Figure [7] The transient nature of the magnetic 
field, and the lack of obvious large-scale structure 
at these contour levels, as compared to the large- 
scale dynamo (Figure[5]), is clear. Nonetheless, the 
mean saturated magnetic field is non-zero; its en- 
ergy is about 4% of the mean energy of the mag- 



netic field in saturation, when averaged over 2.9 
resistive times. Contours of the mean saturated 
magnetic field are presented in Figure [H The field 
is somewhat complex, dominated by to = 0,1,2 
modes, and is weak at the center of sphere. The 
magnetic field at the sphere's surface peaks at the 
equator, though it is composed of to = 1 and 
to = 2 components, in contrast to the large-scale 
dynamo, which only consists of m = 1 modes. 

The growth rate of the magnetic field is much 
larger in this case than in the laminar cases ex- 
amined above. This may be true because the sys- 
tem is much higher above the critical magnetic 
Reynolds number for this set of parameters, as 
compared to the laminar case, or because a differ- 
ent dynamo mechanism is at work, perhaps ampli- 
fication at small scales. If the growth rate of the 
magnetic field were being domi nated by turbulen t 
action at small scales, following Batchelor ( 1950l ). 
we would expect the magnetic energy to grow 
exponentially in time as |B| 2 ~ exp((ei?e) 1 / 2 t), 
where e is the power dissipation per unit mass. 
For this simulation e ~ 16, so for the small-scale- 
dynamo mechanism we would expect the growth 
rate to be ~ 180, not the ~ 40 which is actually 
measured. The critical magnetic Reynolds number 
for this simulation has not yet been found, though 
we have determined that it is somewhere in the 
range Rm < 700. Based on these observations, it 
seems unlikely that the action of small-scale ed- 
dies is the dominant contribution to the growth of 
magnetic energy, but rather since Rm = 2000 the 
rapid growth rate is explained by the high mag- 
netic Reynolds number. 

To gain insight into the length scales important 
to this dynamo, we examine the average angle- 
integrated energy spectra of the velocity and mag- 
netic fields, E(k) and M(k) respectively, where 
M(k) = fc 2 /(|B(k)| 2 )dfi fe , during the growth 
phase of the simulation, normalized to their re- 
spective total energies. These are plotted in Fig- 
ure [9] (the details of how B(k) is calculated are 
given in Appendix |B| . The spectra have a num- 
ber of interesting features. First, the velocity 
field spectrum is very noisy at higher values of 
k. This phenomenon is likely caused by the trans- 
form to /c-space, which involves a radial integra- 
tion across the edge of the sphere, where the ve- 
locity field quickly goes to zero (except for the 
toroidal boundary condition). Since the magnetic 
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0.5 1.0 
Time [t„] 

Fig. 6. — Kinetic and magnetic energies versus 
time for the system forced by the boundary con- 
dition given in Figure [31 with Rm — 2000 and 
Pm = 1.0. 




Fig. 7. — Contours of constant magnetic field mag- 
nitude for four different times during the satura- 
tion of the Rm = 2000, Pm = 1.0 turbulent dy- 
namo. The figures are At ~ 0.2r^ apart in time. 



field does not go to zero at the sphere's edge it 
does not display this effect. The velocity field 
spectrum peaks at relatively low wavenumber, as 
is expected since the energy in the simulation is 
being injected at the largest scales. The magnetic 
energy, in contrast, peaks at higher wavenumber 
and drops off rapidly. The general understand- 
ing of small-sca le dynamos, at least for Pm ~ 1 



and Pm » 1 (jSchekochihin et all 120041 ) is that 



the magnetic spectrum during the kinematic phase 
should peak at k^, which for Pm = 1.0 simulations 
should be the same as k v , the knee in the velocity 
field spectrum. This is not the case in this situa- 
tion, again suggesting that this is not a small-scale 
dynamo. Of course it is worth noting that the 
current understanding of the small-scale dynamo 
exists in the context of a number of assumptions, 
most notably those of homogeneous and isotropic 
turbulence and the lack of a mean flow; none of 
these assumptions apply in this situation. 

Though it has not been ruled out definitely, 
based on the growth rate of the magnetic energy, 
the non-zero mean saturated magnetic field, and 
the spectra of the kinematic magnetic field, we 
believe that the dynamo presented here is not a 
small-scale dynamo, but rather is an example of a 
turbulent large-scale dynamo. It is possible that 
this dynamo is related to some of the previously 
studied turbulent dynamos which po ssess a mean 



flow and simi l ar ma gnetic spectra (jPontv et al 



20051: [m 

lnmni l2006h. or the recen tly-discovered 
'shear dynamo' (jYousef et al. I l2008h . though none 
of these studies were done in a finite domain. 

4. Discussion and conclusions 

The simulations of the proposed experiment 
presented here invoke a number of assumptions 
regarding the forcing at the plasma's outer edge. 
Though we believe that these simulations capture 
the important physics of the experiment, the as- 
sumption that the forcing can be modeled using 
a continuous no-slip boundary condition is debat- 
able. Clearly local effects caused by the discrete 
magnetic and electric fields are being ignored with 
this treatment. Future simulations of the experi- 
ment will endeavor to more accurat ely model the 



forcin g by using the NIMROD code (jSovinec et al 



20041) to add two- fluid, compressible, and collision- 
less physics, including anisotropic and long-mean- 
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Fig. 8. — Contours of constant magnetic field mag- 
nitude for the Rm = 2000, Pm = 1.0 turbu- 
lent dynamo's mean saturated magnetic field. The 
scale has been reduced from Figure [7] for clarity. 



1.e-01 




Fig. 9. — Time- aver aged, angle-integrated, kine- 
matic magnetic and kinetic energy spectra, M(k) 
and E{k) respectively, for the Rm — 2000, Pm = 
1.0 dynamo, normalized to their respect energies. 



free-path effects. Such additions may be especially 
important at lower values of magnetic Prandtl 
number, as one m ight expect the d evelopment of 
Stewartson layers (|Stewartsonlll957fl near the dis- 
crete areas of forcing, similar to effects found in 
cylindrical MRI experiments w ith differentially ro - 



tating upper and lower rings ( Burin et al. 2006). 



One might also expect the plasma to transport 
cusp magnetic field into the volume of the sphere, 
a phenomena that has been observ ed with flow- 
ing liquid sodium ( Volk et al. 2006). This effect 



could be important for rotating convection exper- 
iments, especia lly for the st udy of the physics of 
the tachocline (|Mieschll2005h . 

Since it has long been known that many ax- 
isymm etric flows exist which magnetically self- 
excite (|RobertJl97lHGubbinJl97alDudley & James! 



1989) it is expected that there are a multitude of 
boundary conditions which could be programmed 
into the experiment which would result in laminar 
dynamo action, allowing the study of the different 
means by which large-scale dynamos occur. Other 
boundary conditions which might be of interest in- 
clude solar- type boundary conditions, where the 
equator spins much faster than the poles, and gas- 
giant-type boundary conditions, where the surface 
contains many prograde and retrograde jets. Den- 
sity stratification and convection in rapidly rotat- 
ing plasmas are also possible for this experiment, 
and should be studied. 

The laminar dynamo presented in this study 
is large scale, as must be the case when there is 
no energy in the velocity field at small scales. As 
the fluid Reynolds number is increased the sim- 
ulations become turbulent. As described in Sec- 
tion 13.31 the current understanding of the turbu- 
lent dynamo presented here is that it is not a 
small-scale dynamo, but rather a very turbulent 
large-scale dynamo. Because the outer toroidal 
boundary condition can be controlled in time, it 
is conceivable that a turbulent small-scale dynamo 
could be generated in this experiment without the 
presence of a mean flow. This is a topic which 
requires further study. Another question of con- 
siderable interest is: at what point does a large- 
scale dynamo become so turbulent that it ceases 
to have appreciable energy at the large scales and 
becomes small-scale? Does such a transition ex- 
ist, or does the presence of a mean flow prevent 
this from happening? The proposed experiment, 
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with its ability to reach high values of Rm and 
wide range of kinetic Reynolds number, would be 
uniquely positioned to study these questions. 

It is easily possible to reach Pm 3> 1 under lam- 
inar conditio ns, as been shown in man y numerical 
simulations (jSchekochihin et al.l 120041 ). Another 
goal of the proposed experiment should be to at- 
tempt to study the Pm 1 regime under fully 
turbulent conditions. This will be challenging, at 
least for the boundary condition examined here. 
The forcing strategy of the proposed experiment 
does not fill the volume of the fluid, but only af- 
fects the boundary, and thus a high value of Re is 
needed to make the fluid turbulent. (A simulation 
was performed with Rm = 2000 and Re = 1000. 
This was laminar until saturation of the dynamo 
was achieved.) This need for a high value of fluid 
Reynolds number puts an upper limit on the ex- 
perimental value of Pm that can be reached un- 
der turbulent conditions. Those researchers that 
perform simulations of this physical regime suf- 
fer from a similar problem; for the simulators the 
problem is a lack of the resolution needed to reach 
high Pm at high Rm, for the proposed experiment 
it is one of a technical upper limit on the value of 
magnetic Reynolds number that is experimentally 
accessible. 

In summary, we have presented a concept for 
a high-/? plasma experiment which by its nature 
would be ideal for studying the MHD dynamo and 
a variety of astrophysical fluid dynamics phenom- 
ena, many of which have never been studied exper- 
imentally. The experiment's plasma confinement 
is based on a ring-cusp strategy, and momentum 
is injected into the plasma via differential toroidal 
E X B forcing at the plasma's outer edge. With 
the advances in technology of the last few decades 
this experiment should be able to reach param- 
eter regimes inaccessible to previous multidipole 
experiments. We have demonstrated through nu- 
merical simulations that such a control scheme is 
sufficient for generating velocity fields which are 
capable of both laminar and turbulent dynamo ac- 
tion. The combination of velocity-field tunability, 
high plasma /?, and wide range of kinetic and mag- 
netic Reynolds numbers would make this experi- 
ment a viable choice for exploring a multitude of 
astrophysical phenomena. 
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A. Outer Boundary Values 

The toroidal velocity field radial profile boundary values which are used to generate the results presented 
in this paper, and which are used to generate Figure [31 can be found in Table [3] 

It should be noted that the simulation code uses the following normalization for axisymmetric spherical 
harmonics: 

Y e °(6, 6) = y/(2£+l)P°(cos6). (Al) 
This normalization affects the radial profile boundary values. 



B. Spatial Transforms 



The simulation code evolves the velocity and magnetic fields in a spherical harmonic basis. To calculate 
M (k) and E(k) the fields must be transformed into spherical k space. Rather than convert from the spherical 
harmonic basis to physical space, and then transform to k space, the fields are transformed directly from the 
spherical harmonic basis to spherical k space. This Appendix outlines how this is accomplished. 

The native form of the fields is in terms of radial profiles projected onto a spherical harmonic basis, 
assuming that the fields are divergence free: 



B = ^VXVX [S a (r)Y a {9,6)r}+V X [T a {r)Y a (9, <t>)r] , 



(Bl) 



where we follow the convention of Moffatt <|l978h and use a full r vector, as opposed to the convention 
of Bullard fc Gellmanl (1954) who used the f unit vector. The summation over a is over all valid spherical 
harmonic combinations, I and m, starting at I = 1, truncated at some £ max and m max . The radial profiles 
S a (r) and T a (r) are the profiles which characterize the field. Thus defined, the components of the field take 
the form 



(X 

B e (r,6,<P)=J2 



t a {t a + l)S a {r) 



Y a ( 



1 d(rS a (r)) dY a | T a (r)dY a 
r dr 89 sin# 86 



BJr, 



E 



-r«(r) 



dY a 
89 



1 d(rS a (r)) dY a 



r sin ( 



(B2) 
(B3) 
(B4) 



We now desire 



B(k) = W?T* I = J e^°-B(x)r^sin 



(B5) 



where cos 7 = sin 9k sin cos (6 — 6k) + cos 9k cos 9. Calculation of this integral requires several identities. 
We will use the Rayleigh equation, 



ikr cos 7 



J^i" (2n + 1) i„ (fcr)P„ (cos 7), 



(B6) 



n=0 



where j n (kr) is the spherical Bessel function of the first kind, and P„(cos7) is the Legendre polynomial. 
This expansion over n is truncated at f m ax, which gives satisfactory convergence for the cases considered 
here. We also need the addition theorem for spherical harmonics, 



P„(cos7) = 



4tt 



m=0 



2n + 1 



Yr(e,6)Y^(9kAk). 



(B7) 
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We first consider the radial component of this transform. Combining equations IB21 lB5l IB6I and IB7I gives 



B k (k,9 k ,4> k ) 



n Imax m=n „oo 

-^EE ( n + !) Y n(Qk, A) / Jn(kr)S™(r)rdr, (B8) 

V 27T ' ' Jo 

v n— 1 m— 



where orthonormality between the spherical harmonics has been assumed. The 9 k and <p k components of the 
transform are a little more complicated. Let us consider the first term on the right hand side of equation lB3l 
If the spherical harmonics are defined as Y Q (0,<f>) = N a PP" (cos 9)e ima ^, where N a is the normalization 
constant, then the contribution from the first term is given by 

A Cax ^max m, = U „00 divS™ ) 

v 27T „ . , . „ Jo or 



n=0 e a =l m=0 



dPI n (cos 9) 

P™ (cos 9) — S2jL L S in 9d9. (B9) 

o9 



The second term is similar: 



a £max <m»x m=n „Oo 

Be h (k,e k ,cj> k ) 2 rA A = -jL V V ViV n iV Qi n+1 my n m (^,0 fc ) / Jn (kr)Trjr)r 2 drx 

v n— £ a — 1 m— 

/ P™(cos6)P£(cas6)d6. (BIO) 
Jo 

The terms needed to calculate are similar, the only differences being the quantities in the radial integrals. 

It should be observed that the radial integral is evaluated all the way to infinity. For transforms of the 
velocity field, this integral is only evaluated up to r = 1, since that is where the radial profiles go to zero, 
with the exception of the toroidal radial profiles which have non-zero boundary conditions. (This truncation 
of the radial profile at r — 1 is what gives the velocity field spectra their spiky nature at high k. See Figure [9] 
for an example.) The radial integrals over S a (r) are evaluated all the way to infinity, since the poloidal radial 
profiles are non-zero outside the sphere. Since the region r > 1 is current-free, the radial profile is matched 
to the vacuum solution for the magnetic field. The external part of the radial integrals then become 

/OO />OC 
Jn (kr)S™{r)rdr = S™(1) J j„(fcr)r^"dr, (Bll) 

> d(rS m (r)) C°° 

3n(kr) ^ " rdr = -t a S? a {l) J j n (kr)r-*" dr, (B12) 

the former integral being a special case of the latter. These integrals don't change, and so have been tabulated 
for repeated use. The toroidal part of the magnetic field, T a (r), like the velocity field radial profiles, goes to 
zero at r — 1, since we are assuming an insulating outer boundary. Thus, the radial integrals over T a (r) are 
only evaluated up to r = 1. 
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Table 3 
Toroidal Boundary Values 



I 


Fi{ 


5. EDO 


Fig. [Ib) 


Fig. He) 


Fig. Lid) 


1 


0. 


.5773 


0.2273 


-0.0667 


0.0000 


2 


0. 


.0000 


0.0000 


0.0000 


-0.1002 


3 


0. 


.0000 


0.0541 


-0.1234 


0.0000 


4 


0. 


.0000 


0.0000 


0.0000 


-0.0788 


5 


0. 


.0000 


0.0239 


-0.0130 


0.0000 


G 


0. 


.0000 


0.0000 


0.0000 


-0.0110 


7 


0. 


.0000 


0.0124 


0.0279 


0.0000 


8 


0. 


.0000 


0.0000 


0.0000 


0.0117 


9 


0. 


.0000 


0.0067 


0.0000 


0.0000 


11 


0. 


.0000 


0.0036 


0.0000 


0.0000 


13 


0. 


.0000 


0.0018 


0.0000 


0.0000 


15 


0. 


.0000 


0.0008 


0.0000 


0.0000 


17 


0. 


.0000 


0.0002 


0.0000 


0.0000 
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